Effects of biasing on the galaxy power spectrum at large scales 
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In this paper we study the effect of biasing on the power spectrum at large scales. We show 
that even though non-linear biasing does introduce a white noise contribution on large scales, the 
P{k) cx k" behavior of the matter power spectrum on large scales may still be visible and above the 
white noise for about one decade. We show, that the Kaiser biasing scheme which leads to linear 
bias of the correlation function on large scales, also generates a linear bias of the power spectrum 
on rather small scales. This is a consequence of the divergence on small scales of the pure Harrison- 
Zeldovich spectrum. However, biasing becomes k-dependent when we damp the underlying power 
spectrum on small scales. We also discuss the effect of biasing on the baryon acoustic oscillations. 

PACS numbers: 



I. INTRODUCTION 

One of most promising future observations for cosmol- 
ogy is the precise determination of the galaxy power spec- 
trum. So far, the emphasis has been on the anisotropics 
and polarization of the cosmic microwave background 
(CMB) but in the future we also want to determine the 
matter distribution of the Universe with much better pre- 
cision. The advantage of the matter distribution if com- 
pared to the CMB is that while the latter represents only 
a two-dimensional data set, it has been emitted from the 
surface of last scattering, the former is three-dimensional 
and therefore contains, in principle much more informa- 
tion. The disadvantage is that we can only observe galax- 
ies and it is not clear how the observed galaxy distri- 
bution is related to the underlying matter distribution 
which we calculate using cosmological perturbation the- 
ory (see e.g. [l|) and numerical N-body simulations (see 
e.g. i). 

This problem goes under the name of 'biasing' and has 
been studied in many works. Among the first influential 
papers on the topic are ■ They are based on the idea 
that galaxies form at the peaks of the matter distribution 
or above a certain threshold. A first interpretation of the 
'Kaiser biasing model' was that bias is linear at least on 
large scales But data have shown that bias is not 
linear, see e.g. [1,0]) and this is actually also the case 
in the full 'Kaiser biasing model' [1, Subsequently, 
a fitting formula for bias has been proposed [l^l, which 
is often used, but is not derived from physical consid- 
erations and leads to cosmological parameters which are 
not in good agreement with other observations . Most 
modern versions of biasing have replaced the 'peaks' of 
the density distribution by so called 'dark matter halos' 
inside which galaxies are supposed to form [12 17]. Even 
though these models do qualitatively agree with obser- 
vations [l^, they are very model dependent. One has to 
specify a so called halo occupation distribution (HOD). 
Even sim ple models of this have of the order 10 free pa- 
rameters |19j . Furthermore, simulations suggest that the 



HOD probably depends on time [20|, |2]J and on the de- 
tails of the formation process (e.g. whether the halo un- 
der consideration has been ejected from another halo (2^ 
or undergone mergers [l^l). An alternative idea, where 
galaxies are considered as the fundamental objects and 
halos are attached to them is put forward in Ref. [2^ . 

Another modern approach to biasing is perturbation 
theory. It is based on the assumption that the galaxy 
distribution can be expanded in a Taylor series of the 
density fiuctuation at the same point (local bias) [2^ . 
This can be generalized by allowing also a dependence 
on the divergence of the velocity field and on shear [27| . 
The advantage of this approach is that on large scales, 
where fiuctuations are small, one can reduce it to only 
a few free parameters. The disadvantage is that one 
has no a priori information on the amplitude of these 
unknown parameters and one has to fit them together 
with other cosmological parameters. Furthermore, the 
white noise contribution in the power spectrum on large 
scales is also present in this biasing scheme as soon as 
one allows for non-linear contributions. This has already 
been noticed in [2^ . In addition, the results depend 
on whether one uses Lagrangian or Euclidean perturba- 
tion theory [2^, [s^] . A systematic comparision between 
higher order Eulerian and Langrangian perturbation the- 
ory with numerical simulations can be found in L3l| . 

In this paper we do not develop any detailed biasing 
model, but we want to address two main issues which we 
believe are quite model independent. First, we want to 
illustrate the fact that even though the galaxy correlation 
function may differ from the matter correlation function 
by more than a simple multiplicative constant (non-linear 
biasing) only on small scales, this can significantly modify 
the galaxy power spectrum also on very large scales. 

Secondly we want to estimate the effects of biasing on 
intermediate scales which are relevant for baryon acous- 
tic oscillations (BAO's). These are the remains of the 
oscillations in the baryon-photon plasma generated prior 
to recombination which then contribute to the matter 
power spectrum. These acoustic peaks have been mea- 
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sured in detail in the CMB anisotropy spectrum p33| . 
and there is evidence of their presence also in the galaxy 
power spectrum [s^l • Even though there are some doubts 
about the reality of these features [35*1 , there are now new 
claims of the detection of BAO's even above 4-sigma [36| . 
In any case, substantial progress is expected with future 
surveys like DES [13, BOSS [H or the satehite project 
Euclid Hi]. 

In this brief paper, we are not concerned to use the 
most sophisticated or realistic model of biasing, since we 
mainly want to make a point of principle. In subsequent 
work we shall investigate the ideas presented here in de- 
tail, using a more realistic biasing scheme e.g. the one 
studied in Ref. [Hi. 

In the next section we determine the biased correlation 
function and power spectrum for a simple non-linear bi- 
asing model and discuss the main features. In section Hill 
we conclude. 



II. THE BIASED POWER SPECTRUM 

Non-linear Newtonian clustering does not induce an 
inverse cascade, but rather moves power from large to 
smaller scales (direct cascade). This effect flattens the 
density power spectrum on small scales which accord- 
ing to linear perturbation theory should behave like 
P{k) oc In^(fc) to a behavior roughly like P{k) 
found in numerical simulations, see e.g. On large 
scale, inflation predicts that P{k) cx fc", and the spec- 
tral index n = 0.963 ± 0.014 has been measured with the 
WMAP satellite [s^. This is a very special power spec- 
trum with P(0) = 0. Since the power spectrum is the 
Fourier transform of the correlation function, 

P{k) = \[^ r^{r)joikrydr, j,ix) = , (1) 

V TT Jo X 



as in the entire volume. If this is not the case, one has 
to be cautious of the fact that one is subtracting a ficti- 
tious mean density which is not truly well defined. In this 
sense, here by P{k — 0) we actually mean the theoreti- 
cal infinite volume limit limi;_i.o P{k) in which the power 
spectrum of a sample with no correlations (Poissonian) 
tends to a non-zero constant (white noise). 

We illustrate the point that the biased power spec- 
trum acquires such a 'white noise' contribution by using 
the simple biasing scheme which has been proposed by 
Kaiser [3|: be (5(x) the density fluctuations with vari- 
ance (5(x)^) — a'^ = ^(0). We assume that galaxies form 
when the density fluctuation is larger than a threshold v, 
(5(x) > va. The correlation function of galaxies forming 
according to this prescription is given by the correlation 
function of the threshold sets 9i, defined by [1, 0] 



if (5(x) 
else. 



(3) 



{d{x)) = {9{x)^) = Q{i^) gives the fraction of the volume 
in which (5(x) > h'a. Note that this biasing scheme is 
local, 9^{x.) depends only on 5(x), but, as we shall see, it 
is non-linear. By construction, the amplitude of 0^ never 
exceeds 1 and therefore also its correlation function, 

Ur) = {OA^WAy)) < ^ V r=|x-y|. 

This overall normalization is arbitrary and so is the over- 
all normalization of the corresponding power spectrum 
Pi,. We shall therefore not comment on the overall am- 
plitude of the power spectrum in the following. It has no 
physical significance and could e.g. be multiplied with 
(T^ or (va)"^. 

If we assume that 5(x) is a Gaussian field with vanish- 
ing mean, the one point distribution is 



P{6) 



1 
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P(0) = implies 



^{r)r'^dr = . (2) and the 2-point function is given by 



This integral constraint is very unusual and implies that 
the matter distribution represents a 'super- homogeneous' 
system in a precise statistical sense [io, |4l[ . 

For example, this requires that there is absolutely no 
white noise in the system because this would add a con- 
stant to the power spectrum. From Eq. ([2]) it is clear that 
if biasing is not everywhere linear, leading to a galaxy 
correlation function which is not simply (,g{r) = b^{r), 
there is a high chance that the very subtle and non lo- 
cal integral constraint ([2]) will be violated and Pg{0) ^ 0. 
This actually happens in all the non-linear biasing models 
mentioned in the introduction. Of course, if one subtracts 
the mean galaxy density in an observation to obtain J(x) 
one has, by construction (5(k = 0) = and hence also 
P(fc = 0) = (|(5(0)p) = 0. However, one, in principle has 
to test whether this mean density is well defined in the 
sense that it is the same in half the observation volume 



P{.hM = 



exp 



(4) 



where b\ = (5(xi), ^2 = 5ix-2) and r = |xi — X2|. From 
this it is easy to derive [1, [§1 that the correlation function 
for the biased field Oi, is given by 



Cd.e-^'/^j;:^^^dye-y'/ 



[j^dxe-'/^f 



(5) 



where ^(r) = (z^cr^ — £^{r)x) / ~ ('') ■ Clearly, if 
^(r) = 0, ^(r) = V and hence also $,y(r) = 0, both ^ and 
pass through zero at the same scale tq. Furthermore, 
it can be shown that on large scales, where ^(r) <C 
and (r) <C 1 the correlation functions satisfy the linear 
relation ^i/(r) ~ j/^^(r). For this reason, this biasing 
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scheme is often called 'linear bias'. However, as has been 
pointed out in Ref. 9], the biased power spectrum, 



does no longer vanish at the origin, 



(6) 



(7) 



Even though ^^(r) — i'^£,{r) on large scales, we do not 
have Py{k) — v^P{k) for small k. We now want to es- 
timate the amplitude ^^(O) in order to decide whether 
the turn over to the P{k) oc fc" behavior can be seen in 
the galaxy power spectrum or it is completely masked by 
this biasing effect. Note also, that we have not introduced 
any ad hoc white noise component, but just violated the 
super-homogeneity condition ([2]). As mentioned above, 
this happens for all generic non-linear biasing schemes. 



A. Very large scales 

We have calculated the biased power spectrum P^{k) 
from an underlying ACDM matter power spectrum as ap- 
proximated in [421 fo'^ different values of the biasing pa- 
rameter V. This calculation has already been performed 
in Ref. , but there an exponential cutoff on small scales 
has been introduced for convenience. The reason for us- 
ing such a cutoff was to improve the convergence prop- 
erties of the integrals involved in the computation of the 
biased power spectrum given the underlying one. How- 
ever, as we shall see, such a cutoff can hide some effects 
coming from the small scales that do impact on the re- 
sults obtained in Ref. Here we use the full spectrum 
given in Ref. (4^ which does, however neglect baryons 
and makes use of the fitting formula for the transfer func- 
tion obtained in [5|. Also, this power spectrum neglects 
non-linearities which are in principle relevant on small 
scales. The linear CDM transfer function is given by Q: 

nx) = i^^^i±|^ [1 + 0.284. + (1.18^ 

+ {Q.?,mxf + {O.AQQxfY^'^ (8) 



with X = k/keq and fcgg — •\/2(f -t- Zeq)flMHo- The mat- 
ter power spectrum is given in terms of the transfer func- 
tion by P{k) = T^{k)Po{k), with Po{k) cx fc" the pri- 
mordial power spectrum. This fitting formula does not 
contain non-linear correction to the matter power spec- 
trum that would enhance the small scale power and, with 
it, tj^. This would therefore even enhance the difference 
between the result presented here and in Ref. J9(] . 

In Fig.[5]we show the biased power spectrum obtained 
from the underlying ACDM matter power spectrum. One 
clearly sees that, even though the underlying spectrum 
behaves like fc" at large scales (small fc), the biased spec- 
tra tend to a constant since the integral constraint ([5]) 



is violated, i.e., there is no longer an exact cancellation 
between correlations at small scales and anti-correlations 
at large scales. It is interesting to note that the modifica- 
tion of the power spectrum can be split into two effects: a 
universal modification in the shape and a vertical shift (a 
reduction of power in the present case) that depends on 
the biasing parameter v. This becomes apparent when 
we normalize the biased spectra to their maximum value. 
As shown in Fig [2l the normalized biased power spec- 
tra share the same shape irrespectively of the particular 
value of the parameter v. Moreover, we see that the only 
modification when compared to the underlying ACDM 
spectrum is the appearance of the plateau for very large 
scales, spoiling the P{k) cx fc" behaviour. Thus, we can 
conclude that the biased power spectrum can be factor- 
ized as Pu{k) = 6i(i^)Pi(fc), with bi{i') accounting for the 
i^-dependent shift and Pi{k) the universal biased power 
spectrum shown in Fig. [5] that becomes a constant for 
large scales and tends to the underlying power spectrum 
P{k) on small scales. We obtain the following analytical 
expression for the shift function: 



2cr2 



(I-erf(^./^/2))^ 



(9) 



where erf(a::) denotes the error function. In Ref. "9^ it 
is derived that 6i(z^)^(r) ~ S,u{j') in the regime where 
vS_{r) < 1 and ^(r)/cr^ = £,c{r) ^ 1, hence on sufficiently 
large scales. Here we see that bi {v)P{k) is a good approx- 
imation to the biased power spectrum, Pu{k) on small 
scales as well. This approximation is good for k > k^^ 
up to ^c(l/fc) < 1, hence up to about A: < lO^Mpc"^ 

In the right panel of Fig. [21 we compare the normal- 
ized biased power spectra which all collapse on one line 
with the underlying ACDM power spectrum. Contrary 
to the standard belief, the Kaiser biasing scheme leads 
to linear bias of the power spectrum on small scales and 
non-linear bias on large scales. This might seem to con- 
tradict the well-known fact that biasing is non-linear on 
small scales, as commented in the introduction. However, 
this is not the case. The reason why we obtain the linear 
biasing on small scales here is the log-divergence of the 
pure Harrison-Zel'dovich spectrum on small scales which 
leads to an infinite value of the correlation function at 
the origin, i.e., cr^ = C(0) = oo. Notice that 



/•oo 

cr^ = ^(0) = / k^P 

Jo 



{k)dk. 



(10) 



For large scales (small k) we have that P{k) cx fc" so we 
have no divergences in the lower limit, but for small scales 
(large k) we have P{k) cx (log k)'^k'^~'^ so the integral has 
a logarithmic divergence in the upper limit for n — I. 
We are considering a spectral index slightly smaller than 
1, hence ct^ is not infinity, but very large so that £,c{r) — 
^(r)/(T^ is small already for relatively small scales, see 
Fig. [U and the linear approximation ~ 6i(z^)^(r) is 
valid. 



4 



Of course, on very small scales even the underlying 
CDM power spectrum is damped. Typically, the CDM 
damping scale is expected to be very small, below Ipc. 
However, as we shall show below, once we introduce a 
somewhat larger damping scale which might relevant e.g. 
if a warm dark matter component is present or might 
come from Silk damping of baryons, the divergence of the 
HZ spectrum is regularized, cr is substantially reduced, 
and the amplification of the spectrum on small scales is 
no longer linear. 




r Mpc/h 

FIG. 1: The normalized correlation function = C/i^^ is 
shown. Since the spectral index is close to one, cr^ is very 
large and ^c{r) is small already on small scales. 



We now shall consider the case where small scale power 
has been removed from the underlying power spectrum. 
To that end, we simply smooth the underlying power 
spectrum with a Gaussian window function, 

P = P{k, R) ^ \W{kR)\'^P{k) (11) 

where W{kR) — exp(— i?^fc^/2) is the Fourier trans- 
form of the normalized window function W^(x, R) — 
^exp[-xV(2i?2)]. 

This allows us to compare with the previous results 
obtained in Ref @, since the cutoff introduced there is 
equivalent to the damping introduced here. Also, this 
might be relevant for models with warm dark matter 
where the primordial power spectrum is damped on small 
scales. (It is not clear whether this small scale damping 
survives nonlinear effects which are also relevant on small 
scales and which are neglected in our treatment.) Other- 
wise, this may simple be interpreted as an alternative bi- 
asing scheme: apply Kaiser biasing not to the true power 
spectrum but to the one from which the smallest scale 
structure which is responsible for the very high value 
of a is removed. In this case, the shape of the biased 
power spectrum is no longer universal. It is amplified 
non-linearly both, on scales smaller than the smoothing 
scale and on large scales (see Fig. Indeed, for high 
values of v or large smoothing scale, the turnover that 
is present in the underlying power spectrum disappears 



completely. This disappearance happens for smaller val- 
ues of ly as the smoothing scale increases. These results 
are in line with those obtained in Ref. Q since, as afore- 
mentioned, the smoothing applied here is equivalent to 
the exponential cutoff for small scales introduced in that 
work. As explained above, once we introduce smoothing 
over a scale i?, power on scales smaller than R is damped 
and the UV divergence of the pure Harrison-Zel'dovich 
spectrum disappears. Notice that, cr^, and with it 
strongly depend on the smoothing scale. Therefore, al- 
though the Kaiser scheme leads to linear biasing of the 
non damped CDM power spectrum on small scales, once 
we remove the small scale power, the expected non-linear 
biasing on small scales re-appears. 

If we choose a very small damping scale, like e.g. Ipc 
or even Ikpc, the effect of smoothing is not relevant in 
the observational window from 0.0003/iMpc~^ < k < 
10/iMpc"^ As can be seen in Fig. [3J damping is only 
relevant on scales k^R~^. 

It might be surprising at first to see in Fig. |3] that 
the biased power spectra from larger damping scales are 
higher at given biasing parameter v. But this comes sim- 
ply from the following fact: the larger the damping scale 
the smaller becomes the variance a. For fixed v the prob- 
ability of (5(x)/<T > 1/ then increases and with it the cor- 
relation function and the power spectrum P^. Alter- 
natively, this can be understood by the increase of the 
normalized correlation function = ^/c^ on which 
depends monotonically. Remember however, as stressed 
above, this amplitude is not physical but simply a con- 
sequence of our normalization in Eq. (|3]), 0^ < 1. 

In real observations, the galaxy correlation function, 
i.e. the biased correlation function, has to be smoothed 
over some scale R e.g. with a Gaussian window function 
W{r/R) so that 

= PAk, R) = \W{kR)\^PAk). (12) 

In Fig. 21 we show the corresponding results for the 
smoothed biased power spectrum. We see that the linear 
biasing on small scales disappears, but that the large- 
scale plateau remains because for scales larger than the 
smoothing radius the corresponding spectra are unaf- 
fected so that the non-linear biasing on large scales that 
we have found in the previous cases remain here as well. 



B. The baryon acoustic oscillations 

We now take into account baryons and study the ef- 
fects of the considered biasing scheme on the BAO's. For 
that, we shall use the fitting formula given in [isj for 
the transfer function. Note that the effect of introduc- 
ing baryons is not only the appearance of the BAO's, 
but also a suppression on small scales (Silk damping fT] ) . 
Also here, we use the linear transfer function. 

In Fig. [5] we show our results for this case. In the 
plot we use a too small value of Hm = 2f2f, = 0.1 and 
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FIG. 2: In the left panel we show the underlying ACDM power spectrum (red line) and the corresponding biased power spectra 
(green lines) for p — 5, S, 2 and 1 from top to bottom. As mentioned in the text, the amplitude of the biased power spectra 
is unphysical. In the right panel we have normalized the power spectra to their maximum value after which they all collapse 
to Pi{k) as discussed in the main text. We have also indicated the scales corresponding to the horizon scale at equality, keq 
(blue vertical line), and the present horizon scale, ko (black vertical line). 
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FIG. 3: We show the biased power spectrum for a smoothed underlying power spectrum as explained in the text for several 
smoothing scales (indicated in the corresponding panels) and for u =5, 4, 3, 2 and 1 from top to bottom. Note also here that 
the amplitude of the biased power spectra is unphysical. Actually, the value P(0) is (up to a factor v^) the fraction of the 
volume for which £,c{r) — i.{r)/<j'^ is larger than u which becomes large for larger damping scale R since damping reduces the 



variance a 
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FIG. 4: In this plot we show the biased power spectra after small scale damping. The red line corresponds to the underlying 
power spectrum, the red-dashed line is the smoothed underlying power spectrum and the green lines show the biased power 
spectra after applying smoothing for f = 5,4, 3, 2, 1 from top to bottom. On large scales, k < 1/R, the biased spectra are not 
affected and correspond to the ones shown in Fig. [5] 



h — 0.5 in order to make the BAO's well visible. We 
have performed the calculation also for realistic values of 
the cosmological parameters and found the same results: 
the biased power spectrum flattens at the largest scales 
whereas the shape of the small scales of the spectrum, 
relevant for BAO's, remains unaffected and, therefore, 
the BAO's are not altered. Once again, this becomes 
more apparent when we normalize the spectra to their 
maximum value as shown in the right panel of Fig. [5] As 
before, this result implies that the biased power spectra 
can be factorized as Pv{k) = bi{v)Pi{k) with 61(1^) the 
function given in ^ and Pi the universal shape of the 
biased power spectrum (which of course is different from 
the one without baryons). Hence Silk damping is not 
sufficient to alter the linear amplification of the power 
spectrum on small scales. However, we do find a differ- 
ence here with respect to the case without baryons: the 
power spectrum flattens at a larger scale and its decay 
towards the white noise plateau on large scales is even 
steeper than the P{k) oc A;" behaviour of the underly- 
ing power spectrum. This difference can only be due 
to the different shape of the underlying power spectrum 
on small scales, the baryon acoustic oscillation and Silk 



damping. It is, however not clear whether this finding 
survives in a more realistic biasing scheme. 

Also for this case, we have studied the effect of small 
scale damping of the underlying power spectrum with 
a Gaussian window function. The corresponding results 
are shown in Fig. [51 As before, the shape of the damped 
power spectrum is not conserved after biasing and the 
turnover tends to disappear for high values of the bias- 
ing parameter v and for large smoothing scale. It is in- 
teresting to note that, unlike in the non-smoothed case, 
the BAO's are now distorted by biasing and they are 
completely washed out for high values of u. We find in- 
teresting that smoothing over a scale of only l/i~^Mpc 
already has such a significant effect on the power spec- 
trum on large scales. As in the previous section without 
BAO's, we have linear biasing of the power spectrum on 
small scales which becomes non-linear once we introduce 
smoothing. Again, the reason for this is the high value 
of of the non-smoothed power spectrum which is reg- 
ularized by smoothing. 

Finally, in order to compare to realistic observations 
of the galaxy distribution, we have considered the case 
of introducing smoothing on the biased power spectra. 
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FIG. 5: Biasing in the presence of baryons. For this plot we choose Q,i,h^ = 0.025, Q,Mh? = 0.05 and, h = 0.5. With this too 
large ratio of Q.t,/O.M the BAO's are well visible. In the left panel we show the underlying ACDM power spectrum (red dashed 
line) that includes the effects of baryons and the corresponding biased power spectra (green solid lines) for = 5, 4, 3, 2 and 1 
from top to bottom. In the right panel we have normalized the power spectra to their maximum value that clearly shows the 
two effects explained in the main text. 




FIG. 6: The parameters are like for Fig. [5] We show the underlying ACDM power spectrum smoothed with a Gaussian window 
function (red dashed line) and the corresponding biased power spectra (green solid lines) for i/ = 5, 4, 3, 2 and 1 from top to 
bottom for the different smoothing scales indicated in the panels. 
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FIG. 7: In this plot we show the biased power spectra after applying smoothing. The red line corresponds to the underlying 
power spectrum, the red-dashed line is the smoothed underlying power spectrum and the green lines show the biased power 
spectra after applying smoothing for i/ = 5,4:, 3, 2, 1 from top to bottom. 



The corresponding results are shown in Fig. [T] As in the 
previous section, the non-linear biasing on large scales re- 
mains unaltered because, as commented above, the scales 
larger than the smoothing scale are not affected. It is in- 
teresting to note however that on scales comparable to 
the smoothing scale and smaller, the smoothed biased 
power spectra resemble the smoothed underlying power 
spectrum at small scales. 



III. CONCLUSIONS 

In this work we have studied the effects of a simple 
biasing scheme on the matter power spectrum which has 
been considered before in Ref. [9]. However, we have 
used more realistic underlying power spectra, unlike in 
Ref. where an unrealistic cutoff is introduced. We 
have studied both, the effect of such a cutoff and the 
modification of the BAO's due to biasing. 

First, we have used the approximate BBKS power 
spectrum, which neglects baryons and the non-linear con- 
tribution, and we have found that the biased power spec- 
trum is modified by two effects: a distortion at large 



scales where the power spectrum flattens and a vertical 
shift given by Contrary to the standard claim, in the 
power spectrum there is linear bias on small scales, large 
k, kcq < fc<,10^Mpc~^ and a non-linear (/c-dependent) 
modification of the power spectrum on large scales. The 
linear biasing on small scales is due to the UV divergence 
of the n = 1 Harrison-Zel'dovich spectrum, that makes 
1 already at r = 10~^Mpc for the realistic value of 
n = 0.96 and, therefore, the linear approximation is valid 
already above this scale. We have also studied the effect 
of biasing on a power spectrum where small scale power 
is damped as e.g. in scenarios with warm dark matter. In 
that case, the biased power spectrum becomes distorted 
on all scales. The UV-divergence is no longer present and 
the turnover on large scales tends to disappear for high 
values of the biasing parameter or large smoothing scale. 
Already for a smoothing scale of lft.~^Mpc, the turnover 
on large scales nearly disappears for > 2. If the damp- 
ing scale is very small, R Ikpc, smoothing has no effect 
on the power spectrum on the scales considered here. 

In real observations, there is of course also smoothing 
of the galaxy power spectrum due to finite resolution. 
We have introduced this smoothing in the biased power 
spectrum and we have seen that the linear amplification 
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on small scales also disappears and the power spectra 
flatten on large scales. However, in this case the turnover 
can be visible. 

Finally, we have included baryons and studied the pos- 
sible effects on the BAO's. We have found the same fea- 
tures as in the case without baryons: the biased power 
spectrum flattens at large scales, and it is linearly re- 
scaled by the constant factor 61 (i^) at small scales. There 
are no effects on the BAO's. They are not smeared out 
by Kaiser biasing. However, damping of the underly- 
ing power spectrum distorts the biased power spectrum 
at all scales. For a damping scale of i? ~ IMpc/h and a 
bias parameter f ^2, the BAO's remain intact. However, 
they completely disappear for high values of the biasing 
parameter, z/ > 3 and large damping scale, R > 5Mpc/h. 
Again, we have also introduced smoothing of the biased 
power spectra in order to study the effects in real obser- 
vations. We have seen that BAO's are not affected if the 
smoothing scale is small enough. However, we cannot 
definitively decide on the modification of the BAO's by 
biasing. The result depends on the details of the under- 
lying biasing scheme and is not robust. 

Our findings show that, for the biasing scheme consid- 



ered in this work, the P cx fc" slope is in principle observ- 
able for nearly a decade before the 'white noise' contri- 
bution sets in. When we introduce small scale damping 
of the underlying power spectrum the situation becomes 
worse and the slope can become unobservable, even the 
turn over may disappear. Finally, realistic observations 
require smoothing of the biased power spectra so we have 
also studied this and we have seen that the turnover is 
not affected by smoothing on small scales. This result is 
less pessimistic than in Ref. [§], but it has to be checked 
with a more realistic biasing scheme that will be consid- 
ered in a future project. 
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